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Associated  with  each  arc  (i,  j)  of  a  graph  are  two  numbers  the 

"cost"  and  t^  the  ’'time"  per  unit  flow.  In  our  application  the  unit  flow 
is  a  ship  making  one  trip  from  i  to  J  at  a  cost  c..  and  taking  t 

XJ  lj 

hours.  In  another  example,  a  vessel  for  hire  can  make  a  profit  p^  each 
time  it  goes  from  i  to  J  5  eventually  (if  there  are  a  finite  number  of 
ports)  it  must  complete  a  cycle  with  a  total  profit  P  and  a  total  lapsed 
time  T  where  P  is  the  sum  of  the  profits  and  T  is  the  stun  of  the  times  . 
on  -the  arcs  of  the  cycle.  For  a  maximum  rate  ol  return,  the  shipowners  should 
use  that  cycle  which  maximizes  the  ratio  of  p/t  .  Later  we  shall  describe 
a  more  complex  linear  programming  model  which  we  solve  using  a  column 
generation  scheme  (a  variant  of  the  decomposition  principle).  The 
subproblem  turns  out  to  be  one  of  finding  a  cycle  in  a  graph  that  has  the 
minimum  ratio  of  total  cost  to  total  time. 


Consider  the  following  linear  program: 


Find  Min  z  ,  x^  >  0  such  that 

n 

(l)  2  c..  X..  =  z 

i,  3=1 

n 

(2)  S  t 
i,J=l  J 
n 

xij  =  1 

n 

*ij  - 0 

(5)  s  x.. 
1=1  0 

"  s  xik  “  0 
k<=l 

J  e  1, 2, . . . ,  n 
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Theorem  1:  Associated  with  an  extreme  minimizing  solution  to  (l),  (2),  (3)  is 
a  cycle  whose  total  cost  to  time  ratio  is  minimum. 

Proof:  Let  x^  <=  1  if  (i,j)  is  an  arc  of  some  cycle  and  =  0  otherwise. 
Let  E  t^j  x^  *»  T  ,  then  x^j  =  x^/t  satisfies  (2)  and  (3)  ancl  z  =  c/T  is 
the  ratio  of  total  costs,  C  «=  E  E  c^j  x.^  ,  to  total  time  T  ..  Accordingly  we 
can  always  associate  with  a  cycle  one  of  the  solutions  of  (l),  (2),  (3)  • 

Consider  now  the  class  of  minimizing  solutions  to  (l),  (2),  (3).  We  can  now 
see  that  to  an  extreme  minimizing  solution  corresponds  a  simple  cycle.  This 
follows  because  the  flows  x^  >  0  can  be  represented  as  a  sum  of  simple 
circulations .  If  any  of  these  circulations  had  by  itself  a  lower  ratio 
Ec.  .x  /Et..x  than  another  one,  the  solution  could  not  be  optimal.  Indeed 

1J  xj  lj  lj 

an  improved  solution  could  be  obtained  by  building  up  the  circulation  around 
that  cycle  with  the  lowest  ratio  and  decreasing  the  flow  around  the  one  with 
a  higher  ratio.  Nor  could  a  solution  be  extreme  if  there  were  two  simple 
cycles  with  the  same  ratio  because  one  could  represent  such  a  solution  as  a 
convex  combination  of  two  others  by  first  building  up  and  then  building  down 
the  circulation  in  one  of  the  cycles  while  adjusting  the  other  so  (2)  holds. 

THE  SIMPLEX  ALGORITHM 

A  simple  algorithm  for  solving  (l),  (2),  (3)  can  be  derived  from  the 
simplex  method.  A  basis  involves  n  columns  (one  equation  is  redundant). 

The  corresponding  arcs  in  the  graph  must  consist  of  a  tree  and  one  out-of-tree 
arc.  To  see  this  we  note  that  since  a  basis  is  non -singular,  there  must 
be  at  least  one  non-singular  (n-l)x(n-l)  submatrix  formed  by  deleting  the 
row  associated  with  the  time  equation  and  deleting  some  column  of  the  basis. 

Non -singularity  implies  that  the  n-1  arcs  associated  with  the  remaining 
columns  form  a  tree.  The  arc  associated  with  the  variable  of  the  deleted 
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column  together  vith  a  subset  of  the  arcs  of  the  tree  form  a  cycle. 

In  fact,  this  implies  that  every  basic  feasible  solution  to  (2)  and  (3) 
must  consist  of  one  tree  spanning  all  nodes  and  one  simple  cycle  formed  by 
a  subset  of  the  arcs  of  the  tree  augmented  by  one  additional  arc.  This 
additional  arc  completes  the  cycle  making  possible  the  positive  flow  forced 
by  the  time  equation. 


variables  other  than  the  cycle  variables  have  zero  value  in  a  basic 
solution,  (b)  that  each  node  in  the  simple  cycle  has  one  cycle  arc  pointing 
into  and  the  other  away  from  it,  and  (c)  the  values  of  the  cycle  variables 
are  the  same  and  equal  to  l/T  where  T  is  the  *otal  time  around  the  cycle. 

It  is  also  easy  to  compute  the  simplex  multipliers  (prices)  associated 
with  the  basis.  Indeed  if  ve  let  p  be  the  the  multiplier  associated  with 
the  time  equation  (2)  and  let  IIj  be  those  associated  with  the  node 
equations  (3),  then  for  each  arc  (i,«j)  associated  with  a  basic  variable 

nj  “  ni  +  p  tiJ  =  ciJ 

Summing  these  relations  for  all  arcs  (i,j)  £  Cycle  yields 

(5)  p  =  C/T  =  £  cAj  /£  t^  (i,j)  £  cycle  arcs 
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Knowing  p  ,  one  may  arbitrarily  choose  the  value  of  any  one  node  (there 
Is  a  redundant  equation)  and  determine  the  remainder  by 

(6)  IIj  -  +  (c^  -  p  t^)  (i,j)  €  tree  arcs 

by  branching  out  from  the  selected  node  along  arcs  (i,j)  of  the  tree 
until  all  nodes  are  reached. 

To  obtain  an  improved  solution  the  simplex  multipliers  are  used  to 
eliminate  the  basic  variables  from  the  cost  equation*  The  resulting 
coefficients  for  the  non -basic  variables  are 

(7)  8y  -  (c^  -  V 

If  all  6^j  >  0  ,  then  the  value  of  p  given  by  (5)  is  the  minimum  cost 
to  time  ratio  and  the  problem  is  solved. 

If  not,  let 

(8)  6  <  0  for  some  (p,q) 

Ve  now  make  a  special  Inductive  Assumption:  at  each  iteration,  there  is  a 
basic  feasible  solution  consisting  of  a  directed  tree  spanning  out  from  a 
single  node  that  Is  its  root,  augmented  by  one  additional  arc  to  form  one 
simple  cycle.* 


t  If  a  feasible  cycle  exists  in  the  graph  we  can  satisfy  the  inductive 
assumption  by  taking  any  node  of  the  cycle  as  the  root  of  the  tree  and 
spanning  out  from  the  nodes  of  the  cycle  using  forward  arcs  to  other 
nodes,  and  then  iteratively,  repeating  the  process  with  all  nodes  reached. 

If  nodes  still  remain  they  can  be  reached  by  introducing  high  cost  artificial 
arcs  as  required.  If  there  is  no  feasible  cycle  in  the  graph  there  is 
obviously  no  feasible  solution  to  the  problem;  if  such  is  indeed  the  case 
this  fact  will  be  discovered  by  the  algorithm  that  follows . 
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If  (p,q)  is  the  dotted  arc  marked  "IN  (Case  I)"  in  the  figure,  then  it  is 

easy  to  see  that  entering  (p,  q)  into  the  set  of  basic  arcs  does  not  form 

a  new  cycle  and  we  must  drop  out  of  the  basic  set  the  one  which  is  also 

directed  into  q  .  A  new  basic  feasible  set  of  arcs  is  obtained  with  the 

values  of  II.  decreased  to  II.  +  5  for  node  i  =  q  and  all  nodes  i 

i  i  pq  < 

that  are  followers  of  q  in  the  new  tree.  All  other  II  ^  remain 
unchanged.  In  Cases  II  and  III  a  new  cycle  is  formed  and  we  must  drop 
out  of  the  basis  any  one  arc  (r,s)  which  is  in  the  old  cycle  and  not  in 
the  new  cycle. 

Theorem  2:  If  the  inductive  assumption  holds,  and  if  /(p,q)  is  entered 
into  the  basic  set  in  place  of  the  basic  arc  directed  into  q  ,  then  the 
inductive  assumption  holds  after  the  change  except  when  a  new  cycle  is 
formed. 

Theorem  3:  If  a  new  cycle  is  obtained  as  a  result  of  changes  in  the  basic 
set  of  arcs,  its  p*  =  c/t  ratio  is  less  than  the  previous  one. 

Theorem  4;  If  d^j  =  (c^  -  p  t^)  is  used  as  distances  on  arcs 
(i,j)  in  the  graph,  then  any  cycle  in  the  graph  whose  sum  of  distances 
around  the  arcs  of  the  cycle  is  negative  has  a  lower  c/T  ratio. 

The  simplex  method  accordingly  re  duces  down  to  finding  a  negative  cycle 
in  a  graph  when  arc  distances  d^j  _,lven:  Starting  with  11^  =  0  for  some 

node  of  the  cycle,  the  other  II ^  are  simply  the  distances  from  this  origin 
node  along  arcs  of  the  tree  to  node  i  .  The  simplex  iterative  process  is 
seen  to  be  the  standard  one  for  determining  the  shortest  route  from  the 
origin  to  all  others,  terminating  when  either 
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(9) 


for  all  i,J 


du  -  <nj  -  Vi 0 

or  on  some  Iteration  a  negative  cycle  like  In  Cases  H  or  II  is  found*  In 

the  former  case,  It  Is  seen  by  summing  (9)  for  all  (i,  j)  around  any 

given  cycle  that  distance  around  any  cycle  Is  non-negative  so  that  the 

p  used  to  determine  d^  =  c^  -  p  t^j  Is  the  best  ratio.  3h  the  latter 

case,  a  new  cycle  with  a  lower  value  of  p  Is  obtained  and  new  d, , 

/  *0 

values  are  computed. 


ALGORITHM:  In  the  following  algorithm  the  pairs  ( 1,  J )  represent  directed 
arcs  defined  in  the  graph: 

Or  Let  SQ  be  any  starting  cycle.  If  none  available,  set 
p  =  Max  [c.,/t,.]  In  Step  2  below 

•  lj  lj 


2:  Compute  <=  2  c^/Et^ 
3:  Compute  C1J 


1 :  For  k  =  0,1 , .  •  . 

(i,J)  e  ^ 

for  all  1,  J 

4:  Set  1^  *=  0  and  set  predecessor  of  node  1  as  * 

(meaning  none). 

Set  IT^  «=  «»  and  set  predecessor  for  1/^1  as  1.  f  f 
5 :  For  each  1  <=  1, 2, . . . ,  n  form  for 

J  -  1,2, .  •  >,n  J/^i. 

(a)  If  >  0  return  to  (5)  and  continue  scanning  J  for 
fixed  1  and  then  repeat  Increasing  1  to  1  +  1  until 
i  =  n,  J  *=  n  -  1.  If  i^n,  Jon-1  terminate . 
Cycle  S^  is  optimal. 


-f f These  are  devices  to  initiate  the  computation  without  effort  and  to 
construct  the  starting  directed  tree  necessary  to  satisfy  the  inductive 
assumption. 
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(b)  If  <  0  go  to  (6). 

6:  Determine  the  nodes  In  reverse  order  along  the  route 

R  from  the  origin  to  1  ,  by  back  tracing  the  predecessors 
of  1  . 

(a)  If  J  Is  not  a  predecessor  of  1  ,  then  \change 
predecessor  of  J  to  1  and  replace  value  of 

♦  h)  ,  return  to  (5)  with  1  =  1, 

J  =  2. 

(b)  If  J  Is  a  predecessor  of  1  ,  then  let 

be  the  cycle  along  the  route  R  from 

node  j  to  1  and  back  to  j  along  arc 
(i,j).  Return  to  1  Increasing  k  to  k  +  1. 

For  programming  simplicity  the  above  algorithm  does  not  maintain  a 
directed  tree.  If  It  Is  modified  to  do  so,  the  nodes  can  be  priced 
sequentially  along  the  tree  and  the  return  from  step  6  (a)  to  step  (5) 
modified  to  take  advantage  of  this. 
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COMPUTATIONAL  EXPERIENCE 

Set  I 

Problem 

Nodes 

Arcs 

Seconds  on  IBM  1620 
(including  input -output) 

A 

4 

12 

t 

3.60 

B 

4 

12 

2.16 

C 

5 

20 

3.96 

D 

4 

7 

2.52 

E 

6 

13 

6.84 

• 

Set  II 

(Excluding  input -output) 

F 

i 

5 

20 

1.80 

10 

90 

7.92 

H 

15 

210 

14.04 

I 

20 

380 

33.84 

J 

25 

600 

36.00 

K 

30 

870 

103.68 

For  problems  in  Set  H, 

the  t  values  for 

each  arc  were  randomly 

generated  Integers  between  10  and  60.  Similarly,  the  values  were 
randomly  generated  between  20  and  120. 


We  do  not  have  an  upper  bound  on  the  number  of  operations  except  the 
kind  that  one  could  derive  from  a  standard  proof  of  the  simplex  algorithm. 
In  another  paper  where  a  variant  of  the  scanning  procedure  given  here  is 
used  an  upper  bound  of  (Nodes  +  3)  (Arcs)  additions -comparison  operations 
is  given  for  finding  a  negative  cycle,  [1]. 
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Application  to  a  Ship  Routing  Problem 


Amounts  b^  are  required  to  be  shipped  from  ports  1  to  ports  j  . 
There  are  n  ports  (nodes).  The  shipping  can  either  be  done  by  charter 
at  a  cost  .v^  per  unit  shipped  or  by  using  one  of  a  fleet  of  m  vessels 
under  the  control  of  the  shipping  company.  If  vessel  k  is  used^the  amount 
that  It  carries  between  (i,J)  depends  on  the  kind  of  ship  and  on  the 
pattern  of  ports  forming  a  cycle  g  that  is  assigned  to  the  ship.  We 
denote  this  by  .  Thus  If  arc  (i,j)  Is  not  part  of  cycle  g  ,  then 

a^  »  0  and  If  it  is  its  value  is  the  capacity  of  the  vessel.  trr 

Material  balance  equations :  For  1,  J  =  1,  2, . . . ,  n 

(10) 


n 


yij  +  5  =  1 


x 


kg  c  biJ 


g  Kg 

where  y^  Is  amount  chartered  and  is  the  number  of  times  that  ship 

k  is  employed  in  the  g-th  type  cycle.  We  allow  x^g  to  have  fractional 
values  which  we  Interpret  as  rate  of  use  of  the  ship  in  some  given  period  of 
time. 


Vessel  hour  constraints: 


(ll)  E  +  s^  =  h^  k  e  l,2,...,m. 

g 

is  the  total  hours  available  on  the  k-th  vessel,  is  the  unused 

hours  of  the  chip,  t,  is  the  time  to  complete  one  cycle  of  type  g. 

Kg 

Objective  to  be  minimized: 


where  h^ 


(12) 


z 


i ft  Dependence  on  (i,j)  is  possible  if  type  of  cargo  on  route  (i,j)  is 
different  from  that  on  other  arcs.  In  case  of  airplanes,  capacity  depends 
on  distance. 
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Here  we  are  assuming  that  the  cost  to  operate  vessel  is  per  unit  time 
used,  hence  there  Is  a  savings  of  c^  per  hour  not  used. 

In  an  ore  shipment  application  which  we  were  Interested  In  there  were 
too  many  possible  cycles  to  explicitly  list  all  the  coefficients  of  the 


problem.  Accordingly  we  decided  to  generate  the  column  of  coefficients 

* 

as  needed.  Using  y^  and  sk  as  basic  variables,  one  has  a  starting 
basic  feasible  solution.  We  ndw  assume  we  have  introduced  Into  the  basis 
several  other  columns  and  we  have  a  set  of  simplex  multipliers  p^ 
associated  with  (10)  and  q^  with  (11).  We  wish  to  "price  out"  the 
column  associated  with  x^  and  to  find  that  column  g  for  each  k 
that  prices  out  most  negative.  The  relative  cost  coefficient  of 
becomes 


1.  v  *U  -S  -(1J„  <-*  4  -’■* 

Our  subproblem  becomes  one  of  choosing  that  cycle  in  the  network  of  ship  k 

ilk 

for  which  (15)  is  a  minimum.  Since  a.  J  =  wT.  is  the  ship's  capacity,  if 

Kg  lj 

the  arc  (i,  J )  is  used  in  the  cycle  g  and  zero  otherwise,  the  sum  in  (15) 
is  simply  the  sum  of  the  ship  capacities  on  arcs  (i,j)  weighted  by  p^j 
and  the  times  weighted  by  q^  around  the  cycle  g  .  Note  that  t^  is  the 
sum  of  times  on  arcs  (i,j)  around  the  cycle.  Unfortunately  the  problem 
in  this  form  is  that  of  finding  a  most  negative  cycle  in  a  graph  whose  arc 
distances  are  given.  This  class  of  problems  includes  as  a  special  case 
the  difficult  travelling  salesman  problem. 
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set 


We  got  around  this  difficulty  by  a  change  of  mu 


~  ^kg^^kg*  •®ie  restive  cost  coefficients  for  the  new  problem  become 

<1U>  •  (I  Pij  ^  V 


where  g  denotes  the  (i,,J)  £  cycle  g 

k  ' 

Since  for  fixed  k  is  constant  and  t^  =  L  t^  the  subproblem 
becomes  one  of  finding  that  cycle  g*  that  minimizes  the  ratio 

(1  pij  c«)/(| 


which  fortunately,  as  we  have  seen,  is  a  solveable  probDeml 
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